On the cosmological backreaction for large distance modifications of gravity 
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Every theory that modifies gravity at cosmological distances and that is not already ruled out by 
the Solar system observations must exhibit some nonlinear mechanism that turns off the modification 
close to a compact matter source. Given this nonlinearity, one might expect such a theory to show a 
large gravitational backreaction, i.e. an order one influence of the small scale inhomogeneities on the 
large scale evolution of the Universe. We argue that this is not necessarily the case. If the dominant 
nonlinear terms in the equations obey a shift symmetry, the averaged effect of the nonlinearities can 
be small, although the effect on small scales is large. This happens for DGP (-like) modifications and 
so called f(G) models. For both type of models the shift symmetry is part of the larger "Galilean" 
symmetry. 
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The late time acceleration of the Universe can be in- 
terpreted in two qualitatively different fashions, that are 
both stupefying in their own right. Either we take it as a 
manifestation of some omnipresent dark energy compo- 
nent - dynamical or not- or we view it as an indication 
of the breakdown of Einstein gravity at cosmological dis- 
tances. It is the latter approach that we want to discuss 
in this paper. 

Over the years, in fact already since the early days 
of general relativity, people have been proposing seve- 
ral possible modifications of the gravitational Einstein- 
Hilbert action, ranging from the inclusion of other curva- 
ture invariants or extra fields with non-minimal couplings 
to the introduction of extra dimensions for gravity. The 
corresponding (4D) modified Einstein equations can be 
written generically as: 



(1) 



where the modification appears in the extra tensor H^ v 
that is added to the Einstein tensor at the left hand 
side of the equation. Notice that we are not consider- 
ing modifications to the minimal matter-metric coupling 
in the matter action that gives rise to the right hand 
side. From the uniqueness of spin-2 gravity [l| it fol- 
lows that the modification will necessarily involve one 
or more new degrees of freedom. So if we want to keep 
the equations at most of second order in the derivatives, 
new fields ir will appear in addition to the metric field 
g M „. Actually, in this paper we will restrict ourselves 
to those cases where the dominant contributions to the 
modification arise through only one extra scalar degree of 
freedom, with one corresponding scalar field tt. Variation 
of the action with respect to this new field then results 
in one more equation that we have to consider: 



E w (g,n)=0. 



(2) 



As a first step in establishing the viability of a par- 
ticular model as an alternative for dark energy one then 
looks at the evolution of the scale factor in the homoge- 
nous (presently matter dominated) Friedmann-Lemartre- 
Roberston- Walker (FLRW) universe to see if there exist 



self-accelerating solutions. For those models that pro- 
duce reasonable fits to the supernovae data, the transi- 
tion from an ordinary matter dominated FLRW expan- 
sion to the modified evolution will be controlled by an 
extra dimensionfull constant /i, that is of the order of 
the present Hubble parameter Hq. In any case, the mo- 
dified cosmic evolution implies that on the present cosmic 
background the extra degree of freedom is light, with a 
mass m < Hq , giving it an effectively infinite interaction 
range. 

Of course, the acceleration of the scale factor in the ho- 
mogeneous FLRW background by itself, does not permit 
us to distinguish between large distance modifications of 
gravity and dark energy. It is only when we consider the 
inhomogeneous Universe that we actually live in, that the 
distinction makes sense. By definition, for large distance 
modifications of gravity, the extra light degree of free- 
dom will couple to matter with gravitational strength. 
More explicitly, one can look at the first order terms of 
the full equations in an expansion on the cosmic back- 
ground. Working in Newtonian gauge and dropping the 
vector and tensor perturbations that are subdominant for 
the scalar matter density fluctuations in our Universe, we 
can write the metric for the Universe as (for simplicity 
we will consider a flat background): 



ds 2 



-dt 2 (l + 2i>) + a(t) 2 (l + 2<P){dx 2 +dy 2 + dz 2 ) . (3) 



The rotational symmetry of the FLRW background then 
dictates the following general form for the scalar compo- 
nents of Eqs. |T]) and @ at linear order: 
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where the Oi are linear combinations of <f>, ip and 7r. The 
perturbations are sourced by the matter fluctuations with 
relative density contrast 6 = Ap/p (p is the average mat- 
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ter density). We focus on the non-relativistic, subhori- 
zon density fluctuations, for which the (effective) pressure 
terms can be dropped. Accordingly, we have dropped the 
subdominant terms on the left hand side. For instance, 
we do not write possible mass terms, since as we said, for 
those models that modify the cosmic expansion the mass 
on the present background will be of order Hq. This ex- 
cludes the so called chameleon models from our general 
treatment. For these models the mass of the extra scalar 
field is at least of the order 10 3 Hq, which ties to the fact 
that the late time acceleration is simply driven by the 
cosmological constant term in the scalar field potential 
[3]. Our treatment does however obviously include Ein- 
stein gravity, for which the Ui are zero. In that case the 
first equation is just the Poisson equation for Newton's 
potential <p in comoving coordinates, while the second 
equation results in ip « <p. The parameterization H]) also 
includes DGP (-like), f(R) and /(G) gravity jO, ij. 
It has been shown for instance, that in the subhorizon 
regime the higher dimensional DGP model is effectively 
described by a 4D theory with one extra scalar tt that 
corresponds to the bending mode, which inspired the for- 
mulation of fully 4D DGP-like models [6(. Normalizing 
7r so that 7r ~ <p r*j ijj on the cosmic background, the Oi 
in j4| read: 
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2^±6(1 + ^)tt, (5) 



with the minus sign for the self-accelerating solution 
For f(R) and /(G) gravity, n corresponds respectively to 
the fluctuation of the Ricci scalar and the Gauss-Bonnet 
invariant G |8|. And also for those models the tt fluc- 
tuation mixes up with the metric fluctuations on the 
present cosmic background, resulting in an effective grav- 
itational coupling to matter. This coupling will neces- 
sarily produce a markedly different phenomenology, once 
we bring in the cosmological probes of inhomogeneity 
(CMB anisotropies, large scale structure) in addition to 
the (large red-shift) supernovae observations that probe 
the FLRW background [|. 

Moreover, the gravitational force is not only impor- 
tant on cosmological scales, but also on smaller scales. It 
plays a prominent role in the dynamics of galaxies, stellar 
and planetary systems and as such there are many addi- 
tional observations that can be used to constrain possible 
modifications of gravity . In this context, the tightest 
constraints come from those systems localized on stellar 
system scales or smaller, where we can safely neglect any 
dark components to the energy-momentum tensor. And 
according to the observations of binary pulsars, together 
with the Solar system observations and the table top ex- 
periments on Earth, there is actually very little room for 
any modification of gravity at these distance scales. For 
the PPN parameter 7 « —(p/ip for instance, a careful ob- 
servation of the radio waves emitted between the Cassini 
satellite and Earth, revealed that for the gravitational 



field due to the Sun 7=1 + 0(1O~ 5 ), in agreement with 
the 7=1 result of ordinary gravity In contrast, 

for the modified gravity models that we discussed above, 
the extra scalar polarization in Q on the present cosmic 
background, yields an order one modification for 7, in 
blatant conflict with this constraint. So these models arc 
ruled out if the linearization ([4| holds in the Solar system. 
Conversely, for any viable large distance modification of 
gravity, the Solar system should be in the gravitational 
non-perturbative regime To establish if this is the 
case for a particular model, one has to study the higher 
order corrections to the weak field expansion of H^, E w . 
But let us first briefly review the effect of the nonlinear 
corrections for ordinary gravity. 

For Einstein gravity the dominant higher order correc- 
tions in the quasi-static and subhorizon regime will al- 
ways contain two spatial derivatives. So at second order 
for instance, one has the terms (taking ip as <j>): 
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The local effect of these terms is all too well known. Com- 
paring them with the first order terms of Goo and G\ in 
(|4]), one finds that the linear solutions will break down 
when (p > 1, i.e. in the immediate neighborhood of neu- 
tron stars and black holes. So for Einstein gravity, the 
Newtonianly perturbed FLRW metric ([3]) , with the fields 
solution of the linearized equations (UJ), is nearly every- 
where a good approximation to the actual metric of the 
Universe, also for those systems with highly nonlinear 
density contrasts [13]. Galaxy clusters and individual 
galaxies for instance have 6 ~ 1 — 10, 5 ~ 10 5 respectively, 
while the potential typically has the value <fi ~ 10~ 5 . 

As for the global effect of the terms (|6]), one should 
compare their spatial average with the zero order terms 
G^q , G 1 ^ ~ H 2 that appear in the Friedmann equations. 
This average can be evaluated from the knowledge of the 
density power spectrum. For the modes with wavelength 
L, we have <Pl ~ L 2 H 2 5 from the Poisson equation, so 
that their contribution can be estimated as: 

< 0V 2 >~H 4 L 2 \6 L \ 2 (~4nG N <0 L A PL >), (7) 

with |<5l| 2 the dimensionless spectral density. The last 
term here illustrates its physical meaning as the contri- 
bution to the Friedmann equations of the average gravi- 
tational potential energy of the inhomogeneities. And for 
our present Universe this contribution is of relative order 
1CU 5 which justifies the FLRW framework for cosmology, 
that treats the Universe as approximately homogeneous 
at all scales [l3j |. 

Before looking at the particular models just yet, let 
us now give an intuitive argument, complementing other 
general arguments @, [3], on the type of higher or- 
der terms we would generically expect for cosmologically 
large distance modifications of gravity. Different from 
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Einstein gravity, the modified equations will be nonlin- 
ear functions of curvature. So, in addition to terms of 
the type © that arise in the weak field expansion of cur- 
vature, we should also expect terms that are due to the 
expansion of these nonlinear functions. In particular, on 
cosmic backgrounds with curvature larger or of the same 
order as the characteristic scale of the model, H 2 > pi 2 , 
this expansion will break down for local curvatures that 
are larger than the background curvature. Estimating 
the local curvature due to inhomogeneity as H 2 5, this 
means that the weak field expansion should break down 
for all systems with nonlinear density contrasts. This is 
indeed what happens for DGP and f(G) gravity 15, 
It does not happen for f(R) gravity, the reason being that 
for those models the relevant curvature (R) due to an in- 



homogeneity is locked to its background value [17|. So 
for f(R) models the linearization ([4| holds in the same 
situations as for ordinary gravity, and the contribution 
of the extra light scalar in the Solar system rules these 
models out. 

As we said, for DGP and f(G) gravity, the range of 
validity of the linearization will be different than for Ein- 
stein gravity. From the argument above, we can write 
down the general form of the extra nonlinear terms as 
(taking <f> ~ ip ~ 7r and not being precise about the exact 
derivative structure): 
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Indeed, comparing these terms with the first order terms 
(J4j) , we see that the linearization breaks down for V 2 (f> ~ 
H 2 6 > H 2 . The local effect of these terms for an isolated 
compact object of mass M can be characterized by its so 
called Vainshtein radius R v = (G N M / 'iJ 2 ) 1 / 3 . Within a 
'halo' of radius ~ Ry that is surrounding the object, the 
linear solution for the gravitational field will break down 
[lit Hit ]. Of course, this by itself is not enough to en- 
sure the screening of the extra scalar at short distances. 
However, for DGP gravity, full nonlinear solutions have 
been obtained that match the cosmic background at in- 
finity and for which Einstein gravity is recovered at short 
distances from the object [15j . For /(G) gravity one 
can also argue the recovery of Einstein gravity at short 
distances, although the precise matching to the cosmic 
background at infinity is still an open issue [l8j |. 

The global effect of the new nonlinear terms in our in- 
homogeneous Universe has been largely left unexplored 
so far, and in fact for all studies of the cosmology for these 
models it was (often implicitly) assumed that, just as for 
Einstein gravity, the small scale structure plays a negli- 
gible role. However, both for DGP and f(G) gravity, the 
present Universe consists of compact objects and struc- 
tures separated by distances that are smaller or of the 
same order as their Vainshtein radius. So the lineariza- 
tion will in fact never hold locally and one would not 
a priori expect our Universe to expand in the same fash- 



ion as the corresponding homogeneous matter dominated 
FLRW universe. As a first step in addressing this issue, 
we will only consider those density fluctuations which are 
still in the linear regime, or on the brink of going non- 
linear S < 1. The advantage being that we can still use 
the linear solution that relates the fields to the den- 
sities. In the same way as we did for ordinary gravity, 
we can then estimate the contribution of the modes with 
wavelength L to the average of the nonlinear terms as: 



< V 2 0V 2 > 
H 2 



H 2 \6 L 



(9) 



So we find that already mildly nonlinear density fluctua- 
tions Sl — ► 1 seem to induce an order one backreaction. 
This corroborates our expectations since for those den- 
sity fluctuations the corresponding 'gravitational energy 
density' in ([5]) becomes of the same order as their bare 
density. However, so far, we have not been explicit on the 
precise form of the nonlinear terms. In fact, it turns out 
that both for DGP and f(G) gravity the contributions 
to the averaged equations of terms of the form ([5]) cancel 
out; or that the full average gravitational energy is small, 
of the same order as for Einstein gravity. As we will now 
show explicitly, the reason behind this lies in the shift 
symmetry for the relevant terms in these models. 

For DGP gravity the higher order terms have been 
widely discussed in the literature (see Q and citations 
therein). They are generated by the term in the La- 
grangian (dropping the overall M 2 factor and with our 
normalization for 7r): 



(10) 



Besides having the obvious shift symmetry on 7r, this 
term is also invariant (upon partial integration and drop- 
ping subdominant derivatives of a) under a constant shift 
of the derivatives of ir: 



ir — » 7r + c + c^ ar 



(11) 



This so called "Galilean" invariance is a consistency re- 
quirement, since it prevents the appearance of a new de- 
gree of freedom that has to be a ghost by Ostrograd- 
ski's theorem [19j . Indeed, the symmetry dictates that, 
although the Lagrangian contains four derivatives, the 
equation that follows from it will only contain second or- 
der derivatives, and no new degrees of freedom appear 

However, the normal shift symmetry by itself is suf- 
ficient to ensure a cancellation of the different averaged 
terms. Since by this symmetry the terms in the equation 
En, that follow from the Lagrangian (fT0|) . appear as a 
divergence of the associated Noether current. Explicitly, 
in the quasi-static and subhorizon regime, the relevant 
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terms are: 

= -j^ (V 2 vrV 2 ^ - didjTrdidjn) . (12) 

In the setup ([3]), where the infinite wavelength modes are 
put in the background, the average of a (spatial) diver- 
gence is zero. As one can easily verify, for finite wave- 
length modes, the contributions to the average of the two 
terms on the last line of (fT2|) cancel out. 

The precise form of the higher order terms for f(G) 
gravity has so far not been discussed in the literature. 
We will focus here on those terms that are relevant to 
the central theme of the paper, leaving a more detailed 
discussion for a forthcoming paper [h|. The /(G) models 
are defined by adding some function of the Gauss-Bonnet 
invariant, 

G = R 2 - + R^ pa R^ pa , (13) 

to the Einstein-Hilbert term R of the gravitational La- 
grangian. The choice for this invariant is precisely moti- 
vated by the fact that the corresponding higher derivative 
equations only show one extra degree of freedom in ad- 
dition to the spin-2 graviton :[5|. This can be seen easily 
by eliminating the higher order derivatives through the 
introduction of a Lagrange multiplier A [|| : 

M 2 

C = -g-(R+f'(\)(G-\) + f(\)). (14) 

In this way it is evident that the models are in fact equiva- 
lent to a scalar field with some potential (third and fourth 
term) , no canonical kinetic term and a non- minimal cou- 
pling to G (second term) . Due to the specific total deriva- 
tive structure of G, this non-minimal coupling will gen- 
erate at most second order derivatives (both of A and the 
metric fields) in the full equations @. For those models 
that modify gravity at cosmological distances, the lin- 
ear subhorizon fluctuations of the non-minimal coupling 
term will dominate over the potential. And with a proper 
normalization of the fluctuation it = H 2 f" (Ao)A, where 
Ao is the FLRW background value, the linearized equa- 
tions will then indeed be of the general form Q. This 
non-minimal coupling is also responsible for the domi- 
nant higher order corrections. In the quasi-static sub- 
horizon limit, the breakdown of Q will be due to the 
third order term (again dropping M 2 ): 

4g = /"(Ao)AG^ = -^n (v 2 ^vV - il.nxHl.i),,-) . 

(15) 

We again expect this term to have the full Galilean sym- 
metry, since the equations that derive from it only con- 
tain second derivatives. And indeed, up to boundary 



terms: 

4 tt2 

^-g— 4g = <t> (V 2 V>V 2 7r - didjtpdidjw) 

= 4> (V 2 ttV 2 - didjirdidj^) , (16) 

which shows explicitly that the Lagrangian is invariant 
under a shift of all three fields and their derivatives se- 
parately. The shift symmetry by itself again guarantees 
that the terms, which now appear in all three equations 
Hqq, HI and E v , are divergences that vanish upon aver- 
aging. 

To conclude, our results suggest that a universal shift 
symmetry in the nonlinear Vainshtein screening mech- 
anism, validates the conventional FLRW framework for 
large distance modifications of gravity. However, showing 
that this remains true in a more realistic universe, with 
highly nonlinear density fluctuations i5 > 1, requires go- 
ing beyond our perturbative approach. In this light it 
is interesting to note that the very same shift symmetry 
was used recently to argue the equivalence principle for 
DGP gravity [Hf. 
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